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Abstract 

This paper explores hydrodynamics and hydrostatic of a star in 
post-Newtonian approximation of massive Bran-Dicke gravity. We 
study approximated solution of the field equations upto 0(c -4 ) and 
generalize Euler equation of motion. We then formulate equations 
governing hydrodynamics, stability and instability of the system. Fi¬ 
nally, we discuss spherically symmetric stars for a specific barptropic 
case like dust, cosmic string and domain wall in this scenario. 
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1 Introduction 

The mysteries of the universe always remain interesting issue for the physi¬ 
cists. The combination of two theories, i.e., theory of fluid dynamics and 
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theory of gravity play a major role in the study of the universe. It describes 
universe phenomena in three contexts: (i) dynamics of celestial objects in 
weak-field regimes (like dynamics of solar system), (ii) dynamics of strong- 
field regimes such as Supernovae, formation of black hole, superluminal jets, 
binaries pulsars, neutron stars, X-rays and Gamma-rays burst etc and (iii) 
cosmological issues (origin, different eras and eventual final fate of the uni¬ 
verse) p. 

In this context, the concept of weak-field approximation (Newtonian 
and post-Newtonian (pN) approximations) of relativistic hydrodynamics has 
taken considerable importance in many regards. It provides a platform where 
we can evaluate ranges of deviation and level of consistency between relativis¬ 
tic theory of gravity and Newton’s gravity. The analysis of phenomena in 
strong-field regimes are highly complicated and non-linear, so in order to ob¬ 
tain results of physical interest, weak-field approximation schemes are used 
as an effective tool [2]. Chandrasekhar [3] combined theory of hydrodynam¬ 
ics and general relativity (GR) at pN limits to remove difficulties occurring 
in the analysis of large scale structures. He derived solutions of the field 
equations in pN correction (order of c -4 ) of GR which contain potential and 
superpotential functions to represent masses of gaseous or stellar structure. 
He then used these solutions to derive a set of generalized Newtonian hydro¬ 
dynamics equations and discussed complicated issues of gaseous structures 
(like radial as well as non-radial oscillation of stars, rotating homogenous 
masses and stability of gaseous masses under non-radial as well as radial 
oscillation) in a approximated region in which strong field interactions (grav¬ 
itational radiation) play no role. In this way, he represented a mechanism 
(in weak-field approximation) which can analyze different gaseous or stellar 
systems in GR |Tj. After his work, many researchers [5] explored different 
astrophysical systems in weak-field approximation of GR and obtained many 
interesting results. 

The modified theories of gravity are considered the most fascinated ap¬ 
proaches to resolve mysteries of the present universe ( the dark energy and 
dark matter). Brans-Dicke (BD) gravity [6] is an attractive example of mod¬ 
ified gravity in which gravity is mediated by a massless scalar field 0 and 
curvature. It contains a coupling constant c obd which serves as a tuneable 
parameter and can adjust results according to the requirement. This theory 
has provided convenient solutions of many cosmic issues but unable to ex¬ 
plain “graceful exist” problem of old inflationary according to observational 
surveys. The inflationary model defined by BD gravity is valid for specific 
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ranges of coupling parameter (ujbd < 25) [7] which is not consistent with 
observational limits [8j. Moreover, the BD gravity probes strong field test 
(cosmic issues) for negative and low values of c obd HU but satisfies weak field 
test for high and positive values of ubd [9| - Thus, the weak field evaluations 
are not consistent with the results of strong field. 

In order to resolve this issue, a scalar potential function V(4>) (massive 
function) is introduced in BD gravity which leads to massive Brans-Dicke 
(MBD) gravity [10]. This new theory not only solves cosmic issue (old infla¬ 
tion) but also provides a consistency between the results obtained for weak- 
field as well as cosmic scale [12] - There is a large body of literature which 
describes dynamics of the universe in modified gravity [13] • In this context, 
Nutku [H] explained fluid hydrodynamics in this gravity. He represented 
approximated solutions of BD equations in complete pN limits (order of cT 4 ) 
involving potential functions (representing masses of the celestial objects) 
and explored stability of spherically symmetric gaseous system. Olrno [15] 
evaluated complete pN limits of MBD field equations solutions but he con¬ 
verted only lowest-order (order of c -2 ) approximation of solutions in terms of 
potential functions to discuss f(R ) gravity as a special case of scalar-tensor 
gravity. 

In order to provide a suffice platform for the analysis of gaseous systems 
in accelerated expanding universe, we represent hydrodynamics of fluid in 
complete pN limits of MBD gravity. For this purpose, in section 1, we 
obtain complete pN approximation (order of c -4 ) of MBD field solutions in 
terms of potential and superpotential functions. In section 3 we formulate 
equations governing hydrodynamics of the fluid in pN regimes. Section 4, 
evaluates conditions that govern stability and instability. Section 5 explores 
spherically symmetric barotropic stars in MBD gravity. Finally, section 6 
summarizes the results. 


2 Massive Brans-Dicke Gravity in Post- 
Newtonian Limits 

The action of MBD gravity with (k 2 = ^-) [12j is 

S= ^J + ( 1 ) 
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where L m [g^') represents matter action which depends upon metric as well 
as matter field ■0 / . Variation of the above action with respect to g a p and (j) 
provides respective MBD equations 
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Here T a( a is the energy-momentum tensor of matter, T = g al3 T a p and □ is 
the d’Alembertian operator. Equations (J2J) and (|3]) describe the MBD field 
equations and evolution equation for the scalar field, respectively. We take 
matter contribution in the form of perfect fluid which is compatible with pN 
approximation 

7r 

T a p = [pc 2 (l + -^) + p}u Q up - pg a p, (4) 

where p, /wr, p and u M represent matter density, thermodynamics internal 
energy, pressure and four velocity, respectively. 

In order to calculate some approximated MBD solutions in pN regime, 
the following Taylor expansions are assumed [15] 


9ap ~ Vap + hap, 0 ~ 00 + , 

V{<j>) « Vo + p>Vq + (p 2 Vq /2 + .... 


Here g a p is the Minkoski metric (representing non-dynamical background), 
h a p shows the perturbation tensor which describes deviation of g a p from rj, a p, 
the term 0 O indicates asymptotic cosmic function which slowly varies with 
respect to cosmic time t 0 , (p(t,x) represents local deviation of scalar field 
from (f) 0 with superscripts ((2) and (4)) indicating order of approximation 
(c -2 ) as well as (c -4 ) and Vo = V(cp o) shows the value of potential function 
at time to- The pN limits of MBD theory are obtained through the gauge 
condition 
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and the lowest-order (0(2)) pN approximated solutions describing the po¬ 
tential of the compact object are given by 
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where = 1,2,3 and U is the gravitational potential determined by Pois¬ 
son’s equation 
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The term A 5/5 = 7 + shows the cosmological constant, 7 bd is the parame¬ 
terized pN parameter given by 


_ 3 + 2u BD ~ e~ m ° r 
lBD ~ 3 + 2 lu bd + e ~ m ° r ' 

To discuss characteristics of celestial fluid in complete pN limits (0(4)), 
we evaluate complete pN solutions of MBD equation in terms of potential 
and super-potential functions mm- For this purpose, we calculate values of 
goo ~ 0(4), goi 0(3) in terms of potential and super-potential functions. 
Consider the following metric coefficients 


goo — 1 + hoo, goi — ho i, gij — —Sij + h^, (9) 

where 

hoo = h$ +0(4), h 0i = 0(3), Kj = h\f. 

The components of four-velocity and energy-momentum tensor can be found 
by using Eq.Q given in appendix A. Equation (J2J) can also be written as 


77q/3 — - ~g a pT) + -^p[0, Q 0,p] + ^[0, a; p] - 7^P0 + 7^(0)]- (10) 

The (0, i) and (0,0) components of above equation are given by [15] 
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Here the effect of 0 O is considered almost constant, hence the contributions 
due to 0o and 0 O are neglected. 

I 11 order to solve Eq.flTTT) for hgf, we assume that the second and third 
term of the equation can be expressed in terms of the potential functions x 
as well as U t represented by the Poisson equations as 
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where 10 = —4n Gpi\ [3]. Equations dll]) . (1X31) and (HU) give 
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Now we verify that this solution (HU is consistent with gauge condition as 


h] 


0,a 

dU_ 

dt 


1 h a 

2 a ’° 


d 0 <p 


3 + 2ubd 


3 + 2 u bd + e~ m0r 


3 + 2u B d + e~ mor d 


3 + 2ujbd 


U 


3 + 2 ujbd dxi\ 3 + 2uj BD + e~ m ° r 

From Eqs.OU, (HU and Newtonian equation of continuity, we have 

3 + 2ubd 


(16) 


dU 3 + 2 UBD + e~ mor d 
dt + 3 + 2uj bd dxi V 3 + 2u B d + e~ m ° r 


Ui 


-4v tG. 


eff 


dp d , 

m + alU lv ' 


= 0 , 


6 


































which implies 


dU_ 3 + 2 uj B D + e- mor d f 3 + 2 u B d \ ' 

dt + 3 + 2 u B d dxi \3 + 2 u B d + e~ m ° r ) 


0. 


(17) 


Equations (1T6T) and show the consistency of the solution with gauge 
condition. 

Now we evaluate h[^j from Eq. (fl2l) . For this purpose, we consider that 
the left hand side of Eq. (1T2l) can be defined in terms of potential function 0 
and super-potential given by the following Poisson equations 13 m 
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The complete pN corrections of MBD gravity are given in Eqs. (ED - (d) , (I15p 

and ([20p . Notice that all these solutions are consistent with GR in the limits 

( 2 ) 

Vo -C 1, | <£ t~ 1 and ojbd —> oo. 


3 Hydrodynamics of Fluid in Massive Brans- 
Dicke Gravity 

Hydrodynamics of fluid in any theory is governed by three laws given by 
• Law of conservation of mass, 
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• Law of conservation of momentum, 

• Law of conservation of energy. 

In Newtonian theory, these laws are evaluated through equation of continuity 
and Euler equation of motion. In relativistic hydrodynamics, these law are 
determined through relativistic (generalized) form of equation of continuity 
and Euler equation of motion which are obtained from the following identity 
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3.1 Generalized Equation of Continuity and Euler 
Equation of motion 

According to pN corrections of MBD theory, we have 
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The resulting Christoffel symbols are given in Appendix A. The time com¬ 
ponent of Eq. (I2T1) gives 
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This equation can replace the equation of continuity of Newtonian gravity in 
MBD hydrodynamics. The spatial components of Eq. (12 1 D become 



which can be written as 
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where the potential functions U a . ia , W* and Z^bd) are described in Appendix 
A. Equations ()22i) and (1231) provide generalized form of equation of continuity 
and Euler equation of Newtonian hydrodynamics which represent equation 
of motion in MBD gravity. 

3.2 Conservation Laws 

Here we discuss hydrodynamics of fluid by evaluating the basic conservation 
laws in pN limits of MBD gravity PE!- 

3.2.1 The Conservation of Mass 

The conservation of mass implies that mass neither created nor destroyed. 
The conserved mass can be calculated by integrating Eq. (l22j) over the volume 
occupied by the MBD fluid. The resulting equation of continuity is given by 
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This equation shows that the mass function expressed in terms of density p' 
remains conserved. 


3.2.2 The Conservation of Momentum 

The conservation equations of momentum help us to discuss dynamics of fluid 
motion. The conservation law of linear momentum states that the rate of 
change of total linear momentum is zero, i.e., the total linear momentum of 
the system remains constant. For conserved linear momentum, we integrate 
Eq. (12 3D over the volume along with boundary condition p = 0 so that the 
conservation equation of total linear momentum is 


d 

dt 
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where L t is the total linear momentum per unit volume given by 
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According to law of conservation of angular momentum in homogenous space, 
the total angular momentum of the system remains constant. This is evalu¬ 
ated by first multiplying the equation of motion by Xj and then subtracted 
it by the same expression with indices i and j interchanged. The resultant 
is then integrated over the volume v' giving the equation of conserved total 
angular momentum as 
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3.2.3 The Conservation of Energy 


The law of conservation of energy states that energy neither created nor 
destroyed, i.e., total energy of the system remains constant. In order to 
evaluate total energy of the system, we contract Eq. ([23]) with v, t and then 
integrate over the volume v' . Consequently, the conservation of energy is 
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4 Stability in Massive Brans-Dicke Gravity 

The hydrostatic equilibrium is the state of fluid in which the pressure gradient 
forces are balanced by all other forces (like gravitational forces) and equa¬ 
tion of motion does not depend upon time (z\ = 0). Thus, the hydrostatic 
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conditions from Eq. (l25|) - (127j) are given by 
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where the value of Z^ BD ^ is given in Appendix A. The term “pfl'Bi?” repre¬ 
sents total gravitational effects of MBD systems in hydrostatic equilibrium. 
Equations (l28j) - ([32j) are stability conditions which describe that in stable 
configuration (hydrostatic equilibrium) the total MBD gravitational force 
balances the pressure gradient force due to matter distribution. The total 
density p' due to matter as well as scalar field distribution becomes constant 
(the fluid apparently becomes incompressible). The total linear as well as 
angular momentum of the system vanishes. 

The above discussion implies that in unstable configuration such as grav¬ 
itational collapse, the force of gravity takes over the pressure force (hydro¬ 
static equilibrium of the system is disturbed). The total density of the system 
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does not remain constant and momentum induces into the system. The con¬ 
ditions of such type of instability in MBD fluid can directly be obtained from 
Eqs. ([28]) - (j32j) as 
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5 Spherically Symmetric Barotropic Stars in 
Massive Brans-Dicke Gravity 

Barotropic is a state of fluid in which total density of a system is a function of 
pressure only. In homogenous and isotropic universe (pN regimes), barotropic 
fluid satisfies the equation of state 


P = wp, (37) 

where w is the equation of state parameter. Different values of w represent 
different cosmic configurations such as w = 0 shows dust, w = 1/3 repre¬ 
sents radiation era, w = —1/3 indicates cosmic string, w = —2/3 expresses 
domain walls and w < —1/3 represents dark energy dominated era. It is well- 
known that the configuration of spherically symmetric stars depends upon 
time and radial components only. Here, we represent dynamics of a spherical 
symmetric barotropic stars in pN approximation of MBD gravity. 
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5.1 Hydrodynamics 


Any spherically symmetric star in MBD gravity has the following hydro dy¬ 
namical quantities 


the conserved mass function is defined in terms of density 
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5.2 Conditions of Stability and Instability 

Any star remains stable as long as its hydrostatic equilibrium remains stable. 
In hydrostatic equilibrium, a spherically symmetric barotropic star satisfies 
the following conditions 

% = P3 ' bd (42) 
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From Eqs. (l37|) . (1421) and (t43|) we get 
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The value of Z” (BD ^ is given in Appendix A. 

Since in the phenomenon of gravitational collapse pressure gradient forces 
are over take by gravitational forces. A barotropic star becomes unstable 
(collapses) whenever 
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5.3 Cosmic Strings 


Cosmic strings are 1-dimensional topological defects that are related to soli- 
tonic solutions of the classical equations for complex scalar field. These types 
of matter satisfy equation of state p = —1/3 p [T6]. They are considered to 
have immense density and are significant source of gravitational waves. The 
hydrodynamics of this kind of fluid in MD gravity are described by Eqs. (l38j) - 
(HD along with defined equation of state. The conditions for instability 
configuration (gravitational collapse) are 
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5.4 Domain Wall 

Domain walls are 2-dimensional topological defects in different scalar fields. 
They are topological solitons which are considered as a resultant of spon¬ 
taneous broken of discrete symmetry. These are perfect fluid types which 
obey equation of state p = —2/3 p [1?]. It is believed that collision of two 
such walls violently emit gravitational waves. The equation of state along 
with Eqs.fl3BD-flHD represent hydrodynamics of domain wall in pN regimes 
under the influence of MBD gravity. The respective conditions for unstable 
configuration are 

p < 

p < _g( — f I 9BD dr ) ^ 

3 
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A BD r 2 9 + 6 

3 + 3 + 2 ubd + e~ m ° r 


L r ^ J r g ^ Jr(f>' 7^ 0. 


5.5 Dust Fluid 


Dust is a state of fluid in which fluid particles are either approximately 
stationary or moving in non-intersecting geodesics producing zero pressure. 
Equation (l38l) - (l4Tl) along with p = 0 describe hydrodynamics of dust in 
spherically symmetric configuration. In this case the condition of stability 
(hydrostatic equilibrium) can be directly obtained from Eq. (1421) as follows 


// n 

P9bd ~ 0 , 


since p ^ 0, this implies g" BD = 0 which yield 



where the value of is given in Appendix A. 

6 Conclusions 

This paper is devoted to investigate hydrodynamics as well as hydrostatic of 
MBD fluid in pN regime. For this purpose, we have used simple perturba¬ 
tion scheme and evaluated complete pN approximation of the held equations 
solution in terms of potential functions of celestial masses. We have approx¬ 
imated g 0 o t° 0(c -4 ), 9oi upto 0(c~ 3 ), gij to 0(c~ 2 ) and perturbed scalar 
held upto 0(c -2 ). We have generalized standard Newtonian equation of con¬ 
tinuity and Euler equation of motion in pN limits of MBD gravity. The 
equations governing hydrodynamics, stability and instability of the fluid in 
weak-held regimes are developed. These developed models are then used to 
discuss stability conditions of a spherically symmetric star in pN approxima¬ 
tion of MBD gravity. In particular, we have discussed conditions for unstable 
configuration of barotropic systems such as dust, cosmic string and domain 
walls. 
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The obtained results provide a range of deviation of MBD theory from 
BD and GR gravities in pN approximation. The governing equations involve 
some extra and generalized potential functions that are not used by BD and 
GR theories, which implies that celestial objects used in MBD theory are 
more massive than those described by GR and BD gravities mm- ft is 

also found that MBD gravity reduces to GR theory by freezing the dynamics 

( 2 ) 

at c jbd -+ oo,Vo <C 1 and | |<C 1. The derived model provides an 

essential (boot-strap) character of modified gravity (MBD gravity) which can 
investigate effects of dark energy on the hydrodynamics behavior of large- 
scale systems or we can explore large-scale systems in accelerated expanding 
universe. 


Appendix A 


The contravariant and covariant components of four-velocity and energy- 
momentum tensor are 


u = 


Uq — 


u = 


+ 0(4), 


i + y? 2+w -^f 

1 + ^ ( y 2 - 2U + + °( 4 )- 

A BDr 2 


1 + ^ + 2U 


-) 


— + 0(4), 
c 


ZLi 


Too 

rpOO 

Tqi 

'j-’Oi 


-*? 

rjnij 


-- + 0 ( 3 ), 

c 


= pc 


= pc 


1 + -{v 2 -2U + 


1 H— ?(v 2 ) + 2U — 


pcvi + 0 ( 1 ), 

1 


= pc 


1 + -R+ + 2 U 


Z^-BD 7,2 

3 

2^-bd 7 ' 2 

3 

2^BD r2 


+ 7r) 


+ VT 


+ 0 ( 2 ), 
+ 0 ( 2 ), 


p. 

+ TT + -) 

P J 


+ 0(3), 


T {j = pv l v J + Sijp + 0( 2), 


pv l iP + pSij + — 
& 


’ / 2 , oTv 2 A bdt 2 p, 

p(v +2 U ---b 7r + -)ViVj 

3 p 
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Zm'BoU s t 


2Appr 2 


SijP +0(4). 


The Christoffel symbols in pN corrections are given by 

o _ _ 1 / 3 + 2c jbd , dU 
00 “ _ c3 + 2uj B d + e- m ^ r ’~dt' 

r ° = 1 ® ( 3 + 2ubd tj) 

0i c 2 dxi 3 + 2u B d + e- mor ’ 


T u = 

v 


9^bdU 


dt 

3 + 2c ubd 


+ 4 + 4 ( 


3 + 2cm 


dx.j 3 + 2c obd + e~ m ° r 


+ (A)(_ 3 + 2uJbd _ u.)] _ dx 1 

y dxi y 3 + 2u B d + e~ rn ° r 3 J dtdxidxj \ 


Id. 3 + 2ubd 
c 2 dxi 3 + 2c obd + e _m ° : 


: U) + 


c 4 I dxi \ 2 


- $ --0 -4 


3 T 2c Oi 


dUj + 1 <9 2 y 


3 + 2cj B o + e _m ° r dt 2 <9t<9x 


- lBD-^rS ij ^2 (—)( 


3 + 2c Obd 


dx/ y 3 + 2u B D + e~ mor 


+ (A)(_ 3 ± 2ubd _ u; 

y dxi 3 +2u BD + e~ m ° r 3 ' 


+ ■ 2 ? [sr k (2lBD + 

d Ann 


Appr 


3“ )<5w + df {2lBD + 


X BD‘r- 


d , A BD r 2 

9+ 27bd + - 3 —• 


= 2 = 


djlogV^g) = 2 d 

dx n c 2 &c n 


2 U 7 BD + 


A BD '+ 


which gives 


r 0 M = -0 = -3 


2 a(7BDf/) ^ 9 [„ rr , 

^- Zi -d^ 2XBDU + 


Agp-r 2 

3 


The potential functions U a; i a , W*(x) and Z^bd) expressed in generalized 
Euler equation of motion are defined by 


= O e// / p(x> Q (x') 


, (xj - x')(xj - x-)eh' 
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P 7 

^2 Zi (BD) 


= V ( U i ~ U T,ij) = -Geff j P(X / )^»(X)^(X / ) ( | x _'^3 

- G eff [ p(x') [Vi(x)u Q (x') + V i (x / )v a (x)] <K ' la ^ a ] d * 

+ 3G eff [ p(x') [v«(x)u / 3 (x , )(a; a - x' a ){xfj - x' p )\ 

JV I X X I 

PbdOUbd = \ ( e \ 

c 2 dxi P [ V 3 + 2 u BD + e~ m ° r J 

A BD r 2 2 pv 2 (dU d e~ mor 

3 c 2 \dxi dxi 3 + 2u BD + e~ m ° r 

d(A B pU) dA BD r 2 \ p_ ( dU 


+ - / | 
dxi 


+ jl( 2 ™ 

pc 2 \ dxi 


OXi OXi J pc* \ OXi 

O S e-™ d(A BD U) 1 0A bd C 

dxi 3 + 2c obd dxi 2 dx t 

(tt _ A ^' r2 ^\ dA spr 2 ' 

V 6 ) dxi ' 

The values of Z[ {RD)1 Z'( BD) and Z'" {BD) are given by 


P 7 / _ p'bd 9Ubd _ 

2c 2 i[BD) ~ c 2 dxi ~ P 
p f dU 

+ -4 2 — + 2 - 


P 7 n _ p'bd 9Ubd _ 

2c 2 r{BD) c 2 dxi P 

w ( dU „ d 
+ ^( 2 Ty+ + 2 T: 


ryW 

Zj r(BD ) 


+ *» = ; Li i e ~ mar + A BD r 2 

c 2 dx, t c 2 \ 3 + 2c xbd + e _m ° r / 3c 2 

P_ ( 2 — + 2 - ( — U) + 2+C2+1 

pc 2 \ dxi dxt 3 + 2ubd + e~ m ° r dxi 

I dA BD r 2 \ / A BD r r 2 \ dA BD r 2 ' 

2 dxi J \ 6 / dxi 

P'bd 9Ubd = _i ( e~™ r ,A t^B D r 2 

c 2 dxi c 2 \ 3 + 2uj B d + e~ m ° r J 3c 2 

w ( dU d e -m 0 r di^BpU) 

c 2 \ dr dr 3 + 2u B d + e~ m ° r dr 

ldA BD r 2 \ ( A BD r 2 \ dA BD r 2 

2 dr J + V 6 / dr J ’ 

f 2 / , e" mor r A , A BD r 2 


= L2('i +_++_A + — 

c 2 \ 3 + 2uj B d + e m ° r J 3c 2 

+ 4 _ A BB r 2 \ dA BD r 2 

\ C\ I xl <Y* 
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